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Abstract 



^ ■ We examine the feasibility of generating continuous-variable multipartite entanglement in an 

. intra-cavity quadruply concurrent downconversion scheme that has been proposed for the gen- 

m 

i eration of cluster states by Menicucci et al. [Physical Review Letters 101, 130501 (2008)]. By 

^) ' calculating optimized versions of the van Loock-Furusawa correlations we demonstrate genuine 



quadripartite entanglement and investigate the degree of entanglement present. Above the oscil- 
lation threshold the basic cluster state geometry under consideration suffers from phase diffusion. 



H ' We alleviate this problem by incorporating a small injected signal into our analysis. Finally, we 

investigate squeezed joint operators. While the squeezed joint operators approach zero in the un- 
depleted regime, we find that this is not the case when we consider the full interaction Hamiltonian 
and the presence of a cavity. In fact, we find that the decay of these operators is minimal in a 
cavity, and even depletion alone inhibits cluster state formation. 

PACS numbers: 42.50.Dv,42.65.Yj,03.65.Ud,03.67.Mn,03.67.Lx 
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I. INTRODUCTION 



Cluster states are a class of ffraph states (l| which are of central importance as a re- 

n 

source state for use in one-way, or measurement-based, quantum computing [2|. This type 
of quantum computation, proposed in 2001 by Raussendorf and Briegel [sHsl, differs signifi- 
cantly from the traditional circuit model of quantum commuting in which unitary evolution 
is achieved via a sequence of operations on single qubits 



7|]. In order to realize one-way 



quantum computing, a cluster state is generated and then a sequence of measurements are 



performed on this highly entangled multipartite state 



8|- 



Most quantum computing proposals are based on qubits. Experiments have also been 
performed with qubits. In particular, Grover's algorithm has been implemented using an 
optical one-way quantum computer |9|. This approach relied on generating a four-qubit 
cluster state using a number of independent optical parametric oscillators (OPOs) and beam 
splitters, and then perforrning measurements on this state. 

Lloyd and Braunstein 10| were the first to highlight the potential use of continuous- 
variables (CV) in quantum information. Multipartite CV entanglement has also been ex- 



tensively studied lll-ll4l|. Since then, with the development of one-way quantum computing, 



the notion of using CV cluster states as a potential resource has arisen. Proposals specific to 



CV c 

mm 



uster state quantum computing and the generation of CV cluster states are numerous 



20| . They include schemes based on using a combination of single-mode squeezers and 



quantum non-demolition (QND] gates 171 ] or schemes that rely on single- mode squeezers 



and a network of beam splitters 



181 . This field continues to attract interest and recently an- 



other proposal, based on realizing a CV cluster state using only a single QND gate, was put 



forward 



2l| . Experimental efforts to generate CV cluster states have also taken place, with 



the first quadripartite cluster state generated in 2007 by Su et al. 22| . Similar experiments 
have been performed by Yukawa et al. 1231 . 

n 



Recent efforts |20|, 



24| have also focussed on the possibility of using a single OPO as a 



means of generating CV multipartite entanglement, and in turn, a CV cluster state. In this 
scheme, a single optical cavity is pumped by a number of field modes. The different modes of 
the resonator represent the multipartite entangled systems in the scheme. According to 20|, 
the method generates a cluster state with the quadratures of the optical frequency comb of 
the OPO acting as a quantum computer register. In this article we extend this analysis and 



2 



consider the feasibility of such a scheme. 



In 



2J] a single multimode OPA pumped by two field modes is considered in the un- 



depleted pump approximation. A correspondence is shown between the CV multipartite 
entangled output from this scheme and a CV square-cluster state. It is the square-cluster 



OPO scheme proposed in 



20 



24l | that we consider in our work. Qubit graph states analo- 



gous to the CV square-cluster state have been studied extensively [25| . However, apart from 



26| no study of the entanglement properties of the square-cluster OPO scheme has been 



conducted. Specifically, we investigate this concurrent system in order to verify the pres- 
ence of CV quadripartite entanglement and determine whether or not a CV cluster state is 
produced. This builds on our previous analysis of a similar scheme also based on concurrent 



nonlinear it ies 



271], in which we showed that quadripartite entanglement is present for the 



case of an OPO pumped by four field modes. Furthermore, our work here extends tripartite 



schemes proposed in Refs. [28H30| by utilizing quadruply concurrent nonlinearities. 

This paper is organized as follows. Section II provides an overview of the defining re- 
lation for CV cluster states, presents the Hamiltonian for the scheme under consideration 
and describes the connection between this scheme and the generation of a cluster state. 
In addition. Sec. Ill CI describes the van Loock-Furusawa (VLF) criteria which are used as 
a means of quantifying quadripartite entanglement. Section III considers the interaction 
Hamiltonian in the undepleted pump approximation and gives the VLF correlations under 
this approximation. In Sec. IV we present the full equations of motion for the system and 
calculate the VLF correlations without a cavity present, using the positive-P method. We 
also find the field intensities and compare these to the intensities in the undepleted approx- 
imation. Section V provides an overview of the linearized fluctuation analysis used in this 
work to calculate the measurable output fluctuation spectra from the cavity. These output 
spectra are also found in Sec. V and used to demonstrate violation of the optimized van 
Loock-Furusawa criteria and hence, demonstrate quadripartite entanglement. The spectra 
are obtained above and below the oscillation threshold and the steady-state solutions above 
and below the threshold are also found along with an expression for the critical pumping. 
Finally, in Sec. VI we consider whether or not the entangled output beams produced in 
the proposed scheme do in fact constitute a cluster state, by using the defining relation for 
cluster states and the squeezed joint quadrature operators. 
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II. GENERATION OF A SQUARE-CLUSTER STATE FROM AN OPO 



A CV multi-mode entangled state can be classified as a cluster state if the defining relation 



presented in Refs. 



26| is satisfied. To consider this definition, we first define quadrature 



field operators for each mode as, 



Xi = ai + a], Yi = -i{ai - a|) (1) 

such that [Xj, Yi\ = 2i. We also let X and Y represent column vectors of the amplitude and 
phase quadratures, respectively, for each field mode. The definition of a CV cluster state is 
then any Gaussian state whose quadratures satisfy 

Y -AX — ^ 0, (2) 

where A is the adjacency matrix representing the graph of a given CV state and the arrow 
specifies that the condition holds in the limit of infinite (or large) squeezing. When this 
condition is satisfied for a particular A- matrix, the CV state is a cluster state. The adjacency 
matrix can be weighted and represents the couplings between different nodes on the graph 
representing the cluster state. 



A. Physical Description and Hamiltonian 

The system we model in this paper is comprised of an optical cavity containing a x''^'* 
non-linear crystal. The optical cavity is pumped by two field modes to produce four low- 
frequency entangled output modes at frequencies wa, 074, ws, wg- Mode 1 is pumped at a 
particular frequency and polarization such that it produces modes 3 and 6, as well as modes 
4 and 5. Mode 2 is pumped such that it gives rise to modes 5 and 6. A schematic of the 
setup is shown in Fig. [TJ 

The Hamiltonian for this six-mode system is given by 



'Hint ~l~ Hfree ~\~ 'Hpump ~l~ 'Hbathi (3) 

where the interaction Hamiltonian is 

y-int = ih[xiaialal + Xiaialal + X2a244] + h-c, (4) 
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FIG. 1: (Color online) A x^'^'^ crystal inside a pumped Fabry-Perot cavity. Pump lasers drive 
two intracavity modes with frequencies oji and uj2 (represented by circles and squares), which are 
down-converted to four output modes with frequencies (^3, oj^, oj^ and uq. 



with the Xi representing the effective nonlinearities and hi and a\ denoting the bosonic 
annihilation and creation operators, respectively, for the intra-cavity modes at frequencies 
Ui- The pumping Hamiltonian which describes the cavity driving fields, in the appropriate 
rotating frame is 

2 

j=l 

and the cavity damping Hamiltonian is given by 

6 
i=l 

where ej are the classical pumping laser amplitudes for modes 2, and the Fj are the annihi- 
lation operators for bath quanta, representing losses through the cavity mirrors. 



B. The Undepleted Pump Approximation 

Prior to studying the full Hamiltonian in the presence of an optical cavity, it is useful to 
consider the properties of the Hamiltonian within the undepleted pump approximation. This 
approximation assumes that all the high frequency pump modes remain highly populated 
throughout the interaction process, with no depletion taking place. Specifically, with the 



cavity absent we set = Xi(^i(0)) and ^2 = X2(a2(0)) where are positive, real constants. 
Under this approximation, the interaction Hamiltonian can be written as follows. 



(7) 



]=m,n 



where we assume ^ = = 1, 2), j represents all permutations of the low frequency modes 
and Gmn are the components of the matrix. 



G 



(8) 



1 
10 
10 1 
10 10 

Inspecting the form of this G-matrix for the system, it can be seen that it corresponds 
to the graph in Fig. [TJ The four nodes represent the four low frequency modes and the 
lines connecting the nodes represent the nonlinear coupling of the modes. Such a graph, 
representing the Hamiltonian of Eq. ([7]), has been shown to be equivalent to a four node 
square-cluster state 0, 26 1. 



C. Criteria for Quadripartite Entanglement 

To determine whether or not a square-cluster state generated by the single OPO scheme 
presented in Sec.[TTlis fully inseparable, it is possible to construct multipartite entanglement 
witnesses. These are observables that allow one to distinguish multipartite entangled states 
from separable states. In order to detect CV multipartite entanglement we use the set of 



sufficient conditions proposed by van Loock and Furusawa (VLF) 311], which are a gener- 
alization of the conditions for CV bipartite entanglement 32, IsS]. As shown in these 



conditions may be optimized for the verification of genuine quadripartite entanglement. It 
should be noted that other multi-partite entanglement witnesses also exist js^. 

Using the quadrature definitions in Eq. ([1]), the optimized inequalities which must be si- 
multaneously violated by the low frequency modes in order to demonstrate CV quadripartite 
entanglement are given by. 
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V{X^ - Xs) + ¥{93% + % + % + ge%) > 4, 
V{X^ - Xe) + ¥{93% + ^74^ + ^ + %) > 4, 



(9) 
(10) 

(11) 



where V{A) = (A^) — (A)^ denotes the variance and the gi{i = 3,4,5,6) are arbitrary 
real parameters that are used to optimize the violation of these inequalities. Eq. (Q and 
Eq. ffTOl) are minimized with respect to (74^5 and (73 g, respectively. We then solve the resulting 
equations to obtain the optimized expressions. 



93 - 772 ^717 ' (l^J 

^36 ~ ^3 

= vg^, ' 

= vg^, ' 

V^3(V^46 + ^56) -^36(^34 + ^35) ,s 

96 - 772 T7TP ' (l^J 

•^36 ~ ^^3 

where 

V.^ = !^^^j1±^ -{%){%) (16) 

represents the covariances. For the case where i = j the covariance, denoted V^, reduces to 
the usual variance, V{Yi). 



III. THE HEISENBERG EQUATIONS 



Within the undepleted pump approximation we can calculate the VLF criteria from 
Sec. Ill CI in order to verify the presence of multipartite entanglement. The Heisenberg 
equations of motion for the field operators are given by. 
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dt 

^ = 64 + 64, (19) 
^ = 64 + 64, (20) 

and in turn these equations can be written in terms of the quadrature operators, 



dX, 
dt 

dt 
dX4 
"df 
dn 
~dt 
dX, 

dt 

dY^ 
dt 

dXe 

nr 

dYe 
dt 



^iXe, (21) 

-^1^, (22) 

6^5, (23) 

-^1%, (24) 

eiX4 + 6^6, (25) 

-^lY, - (26) 

^1X3 + 6^5, (27) 

-^1% - ^2%. (28) 



We solve these equations to find analytic solutions for the quadrature operators as functions 
of their initial values, and in turn the optimized VLF criteria can be calculated. Not a great 
deal is learnt from the exact form of these rather complicated analytic expressions. Hence, 
we plot the solutions for the optimized VLF criteria for the cases of equal and unequal values 

We first investigate solutions with both equal. That is, we assume that ^ = ^1,2 and plot 
the optimized VLF criteria. The correlations I36, 145 and I56 correspond to Eq. ([9]) - Eq. (fTTj) . 
respectively. Therefore, a value less than four violates the VLF inequalities. In Fig. |2] we 
see that quadripartite entanglement is present since all three inequalities are simultaneously 
violated. 
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FIG. 2: (Color online) Optimized van Loock-Furusawa correlations, I36 and I45 (black solid line) 
and I56 (red dashed-dot line), found by solving the Heisenberg equations of motion in the un- 
depleted pump approximation. Having all three of the correlations drop below 4 is sufficient to 
demonstrate quadripartite entanglement. All quantities depicted here and in subsequent graphs 
are dimensionless. 

We also consider solutions with unequal C,i, for the case where .^2 = 0.5.^i. The VLF 
correlations are shown in Fig. [31 The violation of I36 and I45 in this case is less than in the 
symmetric case shown in Fig. |2l however, we still observe a substantial violation of these 
VLF inequalities. Furthermore, the violation of I56 is greater here than in the symmetric 
case. 
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FIG. 3: (Color online) Optimized van Loock-Furusawa correlations, I36 and I45 (black solid line) 
and I56 (red dashed-dot line) with ^2 = 0.5^i. 
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IV. THE POSITIVE-P EQUATIONS 



We now turn to an analysis of the setup introduced in Sec. Ill Al by considering the full 
interaction Hamiltonian and introducing a pumped, resonant optical cavity to house the 
nonlinear media. The master equation for this system can be derived using a standard 



approach [35| and is given by, 



dp i 
dt~~h 



Hpump + Hint ,p\ + ^ Ti'^i [p] (29) 



i=l 



where ■ji represent the cavity losses at each frequency and X'Jp] = 2aipa\ — a\aip — pa\ai is 



the Lindblad superoperator [35| under the zero-temperature Markov approximation. From 
the master equation it is possible to derive a set of stochastic differential equations (SDEs) 
and then investigate the intra-cavity dynamics. 



We use the positive-P representation 



36 



38 



to perform a full quantum analysis. In 



this approach the master equation, Eq. (!29|) . is mapped onto a Fokker-Planck equation 
(FPE) for the positive-P function to arrive at a set of c-number SDEs. To obtain the SDEs 
it is necessary that the diffusion matrix of the FPE is positive-definite. In the positive-P 
approach this is achieved by defining two independent stochastic fields and and making 
a correspondence between these operators and the mode operators Sj and a|, respectively. 
This approach allows us to perform stochastic calculations of normally-ordered operator 



moments, for example, in the limit of a large number of trajectories 

Therefore, despite being probabilistic, the positive-P method allows for a full quantum 
treatment of the system when a sufficiently large number of trajectories is used. Following 
this approach, the resulting 12 x 12 diffusion matrix is of the form. 



D = I I , (30) 



08x4 d 

where O^xc are null matrices and d is an 8 x 8 non-zero block given by. 
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J 



(31) 



After factorizing D to find the noise terms, we find the set of Ito SDEs in the positive-P 
representation. The evolution equations for the high frequency fields are, 



dai 
~dt 



£2 - X2tt5a6 - 72^2, 



(32) 



as well as the equations found by interchanging and a^. For the low frequency fields we 
find. 



^ = xittittf^ -73«3 + 'y'^y^(^9(t) + ^^io(i)), 
da^ , 



dt 

da^ 



Xi^ioit + X2a2at - 75«5 + 



(r/5(t) - ir]&{t)) 



da, 
~di 



^ + X2a2a^ - 76«6 + \l^^^{%{t) - iVio{t)) 



+ A/^(r/i(t)-^r^2(t)), 



(33) 
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and also the equations found by swapping with and r]i(t) with 77^+2 (^)- The r]i(t) 
are real, independent, Gaussian noise terms which have the correlations ?7j(t) = and 
Vi{t)Vj{t') = ^ij^{t — t'). We assume throughout this article that all the intracavity modes 
are resonant with the cavity and as a result no detuning terms are included. 

An initial insight into the downconversion processes can be garnered by neglecting the 
pump and loss terms in Eqs. ( 132|) and Eqs. (!33|) momentarily and simply looking at the 
dynamics with depletion present. The results of such a positive-P simulation are shown in 
Fig. m and compared to the undepleted pump approximation results for the low frequency 
modes. Specifically, we plot the intensities of the fields where the horizontal axis is a scaled 
interaction time, with Q = x|cn,2(0)| and x = Xii'i = 1)2). As expected, the undepleted 
pump results deviate from the positive-P results as depletion becomes significant. 

2x10' 
1.5 

0.5 


1 2 3 4 5 

FIG. 4: (Color online) Intensities of the high (red dashed-dot and blue solid lines) and low frequency 
modes (dashed green and small dashed purple lines) calculated using the positive-P equations 
without a cavity present. The number of trajectories is 300,000. The parameters used are xi,2 = 
0.01, 01^2(0) = 1 X 10'^ and 03^4^5.6(0) = 0. The black solid lines show the low frequency modes in 
the undepleted pump approximation. 

We also calculate the VLF correlations using the positive-P equations without a cavity 
present. In particular, we use 300,000 trajectories and all other parameters are the same as 
in Fig. m The results compare well with the undepleted pump results shown in Fig. [2] and 
are not visibly different for a large number of trajectories. 
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V. LINEARIZED FLUCTUATION ANALYSIS 



Returning to our analysis of the more experimentally relevant system introduced in 
Sec. Ill A[ we can consider the full quantum dynamics with depletion present and where 
the interaction occurs inside an optical cavity. This type of scheme is currently under ex- 
perimental study by one of us 



24 



26l | at the University of Virginia. 



We undertake a linearized fluctuation analysis [35| to obtain output spectral correlations 
for the cavity from the intracavity spectra. This is achieved by first linearizing the Eqs. (!32|l 
and fl33|) around the classical steady-state solutions. In the usual manner |35l . 



37|, we then 

find a set of evolution equations for the fluctuations. To begin we neglect the noise terms 
in Eq. (15^ so that af ^ a*. We then set = ai + 6ai, where ai is a mean value and 
6ai represents the fluctuations. This gives a set of classical equations for the mean values 
and from these it is possible to obtain steady-state solutions. It also allows one to obtain 
linearized fluctuation equations from which spectral correlations can be obtained. 

We find that an oscillation threshold is present in our symmetric system. Below this 
threshold we solve the set of classical equations for the mean steady-state values. We find 
that the stationary solutions below the threshold value are 



ai = — for z G {1, 2}, 

li 

ai = fori e {3,4,5,6}. (34) 

Returning to the linearized fluctuation analysis, to first order in the fluctuations the 
equations of motion for the fluctuations. Set = [6ai, 6af, 6a2, . . . , Sa^, Sag]^, are given 
by, 

ddct = -ASfxdt + BdW, (35) 

where B is the noise matrix of Eq. ( l33l) with the steady-state values inserted, dW is a vector 



of Wiener increments 
follows. 



38 



and A is the drift matrix with the steady-state values inserted as 
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where 



(36) 
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VO 72/ 



(37) 
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(38) 
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/ 

For the hnearized fluctuation analysis to be valid the fluctuations must remain small 
compared to the mean values and the eigenvalues of the drift matrix A must have no 
negative real part. The eigenvalues are given by, 
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A3 
A4 

As 
Ae 



71, 
72 - 

72 + 

72 + 

72 + 



1 [^2X2 
2 



72 
£2X2 



71 



V 72 / 



72 



2eixi^2 
71 



V 72 / - 



£2X2 



72 



V 7i / V 72 ^ - 



£2X2 

72 



2eixi 

7i 



£2X2 y 
72 ^ - 



(40) 



There are six other eigenvalues but each of these is degenerate with one of the eigenvalues in 
Eq. ( l40l) . From these expressions it is clear that only A3^4 can have a negative real part. For 
our chosen parameters and in the pump range < ei,2 < 100, it is the eigenvalue A3 that has 
a real part which go from positive to negative. This is depicted in Fig. [5] where we observe 
a plateau of stability for a range of pump values and the transition to an unstable region, 
where the negative real part of A3 is plotted. In the latter region the linearized fluctuation 
analysis is not valid. 

In the following we consider the completely symmetric case where the two pumping 
inputs, ei,2, are equal and given by e, the two nonlinearities, xi,2, are equal and denoted x 
and finally we assume all the cavity losses are equal and given by 7 = 7^ for i = 1, . . . , 6. 
We can now use the expressions for the eigenvalues, as well as the analytic expressions for 
the low and high frequency modes below threshold, to find the critical pumping amplitude, 
ec, at which the oscillation threshold is reached. Furthermore, we confirm via a positive-P 
simulation that this is the threshold pumping value for which downconversion begins to 
populate the low frequency modes. For our symmetric system, we find the critical pumping 
amplitude to be 



11 
X 



(41) 



■1 + V5 

where ec = 61.8 for our chosen cavity parameters. 

If the requirement that the real part of the eigenvalues stay positive is satisfied, the fiuc 



tuation equations will describe an Ornstein-Uhlenbeck process 
spectral correlation matrix is. 



37| for which the intracavity 
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FIG. 5: (Color online) The region of stability (the plateau) and the transition to instability for a 
range of pump amplitudes, ei and €2, found by investigating the behavior of the negative real part 
of the eigenvalue A3, denoted Re[Aj^]. 



S(uj) 



(A + iujl)BB' (A 



ZCijll 



1-1 



(42) 



All the correlations required to study the measurable extracavity spectra are contained in this 
intracavity spectral matrix. Equation fH2]) is related to the measurable output fluctuation 
spectra by standard input-output relations for optical cavities 39|] . In particular, the spectral 
variances and covariances have the general form. 



St,x,(^) = 2V7^5x„x,(c^). (43) 

Similar expressions can be derived for the Y quadratures. For brevity, we use 1°"* (ie. any 
of Igg*, I45*, Igg*) to represcut the three output spectral correlations corresponding to the 
optimized VLF correlations, Ijj, of Eq. ^ to Eq. flTT]) . That is, the same inequalities as 
given in Sec. Ill CI in terms of variances also hold when expressed in terms of the output 
spectra. It is these quantities that can be measured in experiments and that we calculate in 
the remainder of this article. 

A. Output spectra below^ threshold 

We now use these steady-state values to calculate the spectra. In Fig. [H] we plot the 
output spectral correlations, 1°"*, as a function of frequency below threshold for a particular 
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pumping rate. The pump rate is chosen to be e = 0.987ec as this gives the best violation of 
the inequahties for our choice of parameters. The correlations and I45* are equal and give 
the maximum violation of the inequalities. The correlation Igg* is also shown. We observe 
that the spectra bifurcates such that no entanglement is present close to zero frequency. For 
large frequencies — 4. This is the uncorrelated limit for the optimized expressions. 

Collectively, however, these three output spectral correlations confirm that quadripartite 
entanglement is present below threshold, since all I^"* drop below four over a range of 
frequencies. 

5 
4.5 
4 

3.5 
I 3 
2.5 
2 
1.5 

-15 -10 -5 5 10 15 

CO (units of Y) 

FIG. 6: (Color online) The spectral correlations below threshold, for the intracavity quadru- 

ply concurrent scheme. The parameter values are x = 0.01, 7 = 1 and e = 0.987ec. This value of 
e gives best violation of the inequalities. The correlation Igg* is given by the red dashed-dot line, 
while 135* and I45* are equal and given by the black solid line. 




B. Output spectra above threshold 

It is not as straightforward to obtain analytic expressions for the low and high frequency 
modes above threshold. As was seen from the graph in Fig. [1], modes 3 and 4 are only 
coupled to 6 and 5, respectively, while 5 and 6 are also coupled with each other. This 
asymmetry results in an undefined overall phase so that, above the oscillation threshold, the 
low frequency modes suffer phase diffusion, degrading the measurable correlations. We avoid 



this by injecting a small signal 



40 



4l| in one of the low frequency modes, thus stabilizing 



the overall phase. In particular, we find the steady-state solutions numerically using the 
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positive-P equations with an injected signal, 63 = 0.5, in the relevant evolution equations. 
The pump is far more intense than the injected signal, with €3 chosen to be approximately 
one per cent of the pump strength. 

It is also necessary to make use of an injected signal when calculating the output spec- 
tra above threshold. This enables us to calculate the spectra for a range of pump values 
above threshold. In Fig. [7] we plot the output spectral correlations above threshold for a 
pumping rate of e = 1.49ec. We see that all three of the VLF inequalities are violated for a 
range of frequencies and therefore quadripartite entanglement is present above the threshold 
condition. As in the below threshold case, we also see the spectra bifurcating close to zero 
frequency and hence not demonstrating entanglement in this region. 

5 
4.5 
4 

3.5 

3 
2.5 

2 
1.5 

-15 -10 -5 5 10 15 

CO (units of Y) 

FIG. 7: (Color online) The spectral correlations above threshold for the intracavity quadruply 
concurrent scheme. The parameter values are x = 0.01, 7 = 1 and e = 1.49ec. The other 
parameters are the same as in the below threshold case. The correlation is repesented by the 
red dashed-dot line. The correlations and I45* are equal and represented by the black solid 
line. 

We also calculate the maximum quadripartite entanglement for a range of pump field 
amplitudes below and above threshold. This result is shown in Fig. [HI where we plot the 
minimum value of the output spectra at any frequency, as a function of e/ec- From this plot 
it is clear that quadripartite entanglement persists below threshold and well above threshold. 
In Fig. M and Fig. [TD] we plot the VLF correlations as a function of frequency and pumping 
rate. In Fig. [HI we crop the peak of the spectra for visualization purposes. The value of this 
correlation would otherwise increase to a maximum value of I36 45 ~ 60. From these surface 
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plots we can also see the onset of the bifurcation. 




0.5 1 1.5 2 

e/e 



FIG. 8: (Color online) Maximum quadripartite entanglement as a function of the ratio of the 
cavity pumping to the pumping threshold. All other cavity parameters are the same as in Fig. [6] 
and Fig. [71 Again, the correlations V^q* and I45* are equal for these parameters and hence cannot 
be differentiated on the graph. Note that for e/ec = 1 the validity of the results is limited as the 
linearized analysis is no longer valid. 




FIG. 9: (Color online) Plot of the spectral correlations I§g* and I45* as a function of the pumping 
rate e/ec and frequency u (in units of 7). All other cavity parameters are the same as in Fig. [HI 
Note that we crop the peak of the spectra for visualization purposes. 
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FIG. 10: (Color online) Plot of the spectral correlation as a function of the pumping rate e/ec 
and frequency oj (in units of 7). All other cavity parameters are the same as in Fig. [HI 

VI. CONFIRMING THE CLUSTER STATE DEFINING RELATION 

Now that we have demonstrated genuine quadripartite entanglement below and above 
threshold, we turn our attention to whether or not the state generated by the scheme is in 
fact a cluster state. To do this we employ the defining relation for cluster states, given by 
Eq. ([2]). In the undepleted pump regime, there exists a connection between this condition 



and the squeezed joint operators for the system [23|, |2J, |26|, |42| . Under certain mode rotations 
and in the limit of infinite (or large) squeezing the operators for the system give rise to the 
cluster state equation of Eq. (E]) for a square-cluster state. We stress that this equivalence has 
only been shown to hold within the undepleted pump approximation. We briefly overview 
the equivalence here. The squeezed joint quadrature operators are simply the eigenvectors of 
the system found by solving the Heisenberg equations of motion for the system represented 
by the Hamiltonian in Eq. ([7]) with the adjacency matrix given by Eq. (|8]). For our scheme 
the joint quadrature operators are. 



01 = (-ClX3 + ClX4-X5 + X6)e-^^^ 

02 = (-C2X3-C2X4 + X5 + X6)e-^l^ 

03 = (C1F3 + C1F4 + ^^5 + %)e-''\ 
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where Ci = (-\/5 — l)/2, C2 = {V5 + l)/2 and r = is the squeezing parameter. The 
common eigenstate of these joint operators is a quadripartite entangled state that tends 
towards a cluster state when r — oo. That is, any squeezing operator combinations that are 
proportional to the squeezing factor e~^' will automatically satisfy the cluster state condition 
of Eq. ([2]) in the limit of infinite (or large) squeezing. We investigate this by example and 
compare these operators to the cluster state defining relation, Y — AX — > 0. We choose 
one possible solution for the adjacency matrix A, given by 



A = - 
2 



1 V5 

V5 1 

1 75 
V5 1 



(45) 



and which corresponds to a weighted square graph CV cluster state. From Eq. (|2]) the 
resulting cluster state equations are. 



•^2 2 
X, V5X, 

(46) 

where the arrow again represents the limit r — )■ oo. In the undepleted regime and in the 
limit of infinite (or large) squeezing these equations are equivalent to Eqs. f l44l) under certain 
mode rotations. Specifically, if we rotate modes 5 and 6 by 7r/2 by substituting X ^ Y and 



Y — 7- —X the equivalence can be seen 2J]. Hence, by simply verifying that the squeezed 
joint quadrature operators approach zero in the limit r — t- oo, we can determine whether or 
not the proposed scheme gives rise to a cluster state. As an aside, we actually calculate the 
variances of the squeezed operators of Eqs. (1441) . and ensure that these approach zero, as 
these are the quantities we have access to in our simulations. 



We first confirm that the squeezed operators approach zero as expected 2J] in the unde- 
pleted case. In Fig. [11] we calculate the squeezed operators from the Heisenberg equations 
(solid lines). We see that Oi = O3 and O2 = O4 and both sets of operators approach zero as 
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r = approaches infinity. Therefore, the cluster state equation is confirmed in the unde- 
pleted regime. In Fig. [TT]we also plot the squeezed operators calculated from the positive-P 
equations in the absence of an optical cavity (dashed lines). We note that the positive-P 
results go to zero but start to increase from zero at longer times. Moreover, as 02,4 ap- 
proaches zero, 3 increases and does not approach zero. For the cluster state relation to 
be satisfied, all squeezing operators must approach zero simultaneously. As seen here this is 
not the case, except at around C,t = 2.5 for the system considered, and hence the required 
cluster state condition is only partially satisfied once depletion is included. Thus, depletion 
alone is enough to inhibit cluster state formation. 



10 

8 
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FIG. 11: (Color online) Squeezed joint operators, Oj, in the undepleted pump approximation 
(black solid lines) and positive-P results without a cavity (colored dashed lines). The number of 
trajectories for the stochastic simulation is 300,000 with x = 0.01 and ai,2(0) = 1 x 10^. 

In Fig. [12] we plot the squeezing operators for the full interaction Hamiltonian in the 
presence of a cavity, based on positive-P simulations. This result is more relevant to the 
scheme that would be realized in the proposed experiment. Unlike the undepleted case, the 
squeezing operators calculated in this case do not approach zero. Instead, the squeezed oper- 
ators plateau at non-zero values in the steady-state. We find that the decay of the squeezed 
joint operators is in fact quite minimal (of order 20%). This indicates that generating a 
cluster state from the output of a single OPO may present a challenge. 
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O 6 



t (units of 7 ) 

FIG. 12: (Color online) Squeezed joint operators, Oj, calculated using the positive-P equations in 
the presence of a cavity. The operators are calculated for the below threshold case using 50,000 
trajectories, with e = 0.6472ec, x = O-d and 7 = 1. 

VII. CONCLUSIONS 



We have examined an experimentally feasible quadruply concurrent intracavity scheme 



as proposed in Refs. [20|, |2J|. We investigated this system as a potential source of CV 



quadripartite entanglement and as a candidate for the generation of a CV square-cluster 
state. We verified the presence of quadripartite entanglement via optimized versions of the 
well-known VLF correlations. The proposed scheme provided a source of bright entangled 
output beams above the critical pumping threshold when an injected signal was incorporated 
into the analysis. Below threshold we also detected quadripartite entanglement with the 
maximum entanglement predicted near threshold. 

We have also calculated the squeezed joint operators to determine if the state produced by 
the proposed scheme is a CV cluster state. Within the undepleted pump approximation, we 
confirmed that the squeezed joint operators approached zero in the limit of large squeezing. 
However, in our analysis of the more experimentally realistic case where depletion is present 
and a cavity is used to house the nonlinear media, we did not observe that the squeezing 
operators approached zero as required by the cluster state defining relation. Furthermore, 
including depletion alone was sufficient to inhibit cluster state formation and once the cavity 
was also included the decay of the squeezing operators was found to be minimal. This leads 
us to conclude that the utility of this system as a source of cluster states depends on the 
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degree to which less than perfect squeezing is acceptable. Overall, solution of the undepleted 
Heisenberg equations of motion can only be used as a general guide to the performance of 
such a system once it is placed inside an optical cavity. Without the cavity, we also find 
that an energy- conserving positive-P simulation of the basic downconversion process shows 
inhibited cluster state formation due to pump depletion. 
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